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ABSTRACT 

Most massive stars are found in the center of dense clusters, and have a companion fraction much 
higher than their lower mass siblings; the massive stars of the Trapezium core in Orion have ~ 1.5 
^sO ' companions each. This high multiplicity could be a consequence of formation via a capture scenario, 

or it could be due to fragmentation of the cores that form the massive stars. During stellar formation 
circumstcllar disks appear to be nearly ubiquitous. Their large radii compared to stellar sizes increase 
the interaction radius significantly, suggesting that disk interactions with neighboring stars could 
assist in capturing binary companions. This mechanism has been studied for stars of approximately 
solar mass and found to be inefficient. In this paper we present simulations of interactions between a 
22 Mq star-disk system and less massive impactors, to study the disk-assisted capture formation of 
, binaries in a regime suited to massive stars. The formation of binaries by capture is found to be much 

■ more efficient for massive capturers. We discuss the effects of a mass dependent velocity dispersion 

and mass segregation on the capture rates, and consider the long term survival of the resultant binaries 
t-H ■ in a dense cluster. 

' Subject headings: circumstellar matter — stars: formation 

m : 
m . 

^ \ 1. INTRODUCTION 

The high multiplicity of massive stars ijMason et alJ Il998t iPreibisch et all 119991: iStahler et all 120001 
\Q ! iCarcia fc Mermilliodl 120011 may provide clues to th eir formation. A useful number to describe the multiplicity is 
O ■ the companion star fraction l|Zinnecker fc Batdl2002T) . csf = (B + 2T + 3Q + ...)/(S + B + T+ ...), with S the number 
of single stars, B the number of binaries, T the number of triples, and Q the number of quadruples. In the Orion 
Nebula cluster, the Trapezium stars have csf = 1.5, high compared to csf ~ 0.5 for low mass stars in the cluster. 
Fragmentation and capture are the two most likely mechanisms by which these binary and higher multiple systems 
form. Capture could occur dynamically, via a three body encounter, or by dissipation of orbital energy through disk 
+3 . or envelope interactions. 

Capture by a disk has been studied as a binary formation mechanism by previous authors. iClarke fc Pringlel l)1991f ) 
performed an analytical approximation to derive an upper limit to the capture rat e for roughly 1.0 Mq star-disk 
systems, concluding that the rate is too low to contribute to the observed multiplicity. iHellerl l)1995|) used a smoothed 
particle hydrodynamics (SPH) code to study encounters between 1.0 M stars, one of which is surrounded by a 100 
AU, 0.1 Mq disk. These conditions, when the envelope is depleted and the prestellar material is found solely in a disk 
around the star, represents the class II phase of prestellar evolution. He examined a range of inclination angles and 
periastra, and from these results determined a binary formation rate. The rates found for this late stage of evolution 
a re too low to implic ate capture as a significant binary formation mechanism. 

Boffi n et alJ l)1998f > performed a similar analysis to Heller, but on a less evolved, more extended disk system. The 
impactor was again 1.0 M , but the primary and disk were both 0.5 Mq, and the disk radius was 1000 AU, similar 
to a class object. The motivation for studying a less evolved system is th e observation th at by the time stars reach 
the class II and III stage of evolution, their binarity has been established l(MathieiJll994fl . Even with the increased 
cross sections of these larger disks, the capture rates in dense regions like the Orion Trapezium center are too low to 
account for the observed multiplicity fractions. 

The early evolution of massive protostars is less well understood than lower mass cases. Disk or envelope capture 
during the formation of a massive star could make a significant contribution to their observed multiplicity. Because 
of its stronger gravitational focusing, a 20 M star with a 2 M disk is a more attractive encounter partner for a 
lower mass cluster sibling. When encounters occur, the interactions will involve a disk with a mass comparable to the 
impactor. The lowest mass impactors will encounter a disk an order of magnitude more massive than themselves. The 
r egime with secondary st ars significantly less massive than the disk has not been explored. 

iMoe ckcl & Ballv ( 2006) performed simulations of encounters between a massive star-disk system and a less massive 
intruder. That work considered encounters at only one radius and impactor mass, focusing on the effects of repeated 
encounters on the disk. In this paper we extend the analysis of encounters with massive protostars, and describe the 
results of simulations of encounters between intruder stars of several different masses and a 20.0 M star surrounded by 
a 2.0 Mq disk, 500 AU in radius. Though the applicability of the standard prestellar classification scheme to massive 
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protostars is unclear, this moderately-extended disk is similar to a later class II stage in the formation of a massive 
star. 

Th ere are several ob servational examples of disks around ma ssive protostars, e .g. NGC 7538 IRS 1 /2 llKraus et alJ 
[200^ . G192. 16-3.82 iShepherd et al.ll200l. Cepheus A HW2 llPatel et al.ll2005^ . IRAS 20126+4104 l|Cesaroni et alJ 
Il999li . and MWC 349 l|White fc Beckerlll985l iTafova et alJl2004l) . A commonality among these sources is that they 
are all young, embedded protostars. By the time a massive star has reached a more evolved stage, for instance in the 
Trapezium, disks are absent. Either during or shortly after the embedded phase, disk destruction is apparently fast, 
via photoevaporation or other mechanisms. The disk lifetime is then approximately the length of the embedded phase, 
~ 10 5 — 10 6 yr. We show that over similar time periods, encounters with cluster siblings can contribute both to the 
formation of binaries and the destruction of the disks around massive stars. 



2. SIMULATIONS 

Th e code used in this study is a modified version of the publicly released SPH/N-body code GADGET-2 ijSnriugel 
2005) . We model the stars as point masses interacting only through gravity, while the gas particles experience 
hydro dynamical forces as well. We have modified the code so that the stars behave as sink particles l)Bate et alJ 
1995), accreting gas that falls within the pre-set, mass dependent accretion radius. The accretion radii are set to be 
small compared to the Bondi-Hoyle radius r g ~ 2GAI+/ (v 2 + c 2 ) 1 ^ 2 , where c s is the sound speed. Using sink particles 
prevents the gas from clumping in unphysical manners about the stars, and prevents very dense gas from dominating 
the integration once it is very close to the stars. We have a lso modified the numerical viscosity slightly, to remain at 
lo w levels except when shoc ks are present ijMonag han 1997). For more detail on the code and the changes made to it, 
see lMoeckel Rr. Ba.llvl l|200^ . 

We set up the disk with a surface density profile 

E(r) = E , (1) 

with the vertical density structure at a given radius 

p(r,z) = p (r)exp 



2H{rf 



where H(r) is a temperature dependent scale-height and po the density at the mid-plane. The smoothing lengths of 
the SPH particles are less than H(r)/2 for all but the innermost regions of the disk. 
The temperature profile is 

f r\- 1/2 

T(r)=Ib(-) , (3) 

with To set depending on the mass of nearest stellar particle. To return to this temperature profile when the disk is 
perturbed, we implement a simple cooling scheme, 

du _ U - Ub ase 
dt t c 

Here r c is the cooling timescale, and is inversely proportional to the Keplerian orbital frequency at the particle's 
radius. The temperature profile of equation © is converted to an internal energy to provide Ubase- The disks are not 
allowed to cool below the base temperat ure profile; in the absence of this , the cooling timescale we use is short enough 
that fragmentation would likely occur ijGammi e3 l200ll lR ice et al.ll2005j) . an effect we are not considering here. The 
temperature and density profiles are appropriate to a disk shortly after its formation ijLin &: Pringldfl9 90K 

We place ~ 3.2 x 10 4 particles pairwise symmetrically around the central star, and normalize the surface density so 
that 2 M Q is contained in the disk out to r<2 = 500 AU. The temperature is normalized so that the disk is Toomre stable 
throughout. We allow the disk to evolve in isolation until it has reached a steady state, at which point we introduce 
the impactor. The initial disk is shown in the first panel of figure ^ To test the resolution of our simulations, we 
ran a retrograde encounter with a 1 M impactor, and periastron radius 0.6 times the disk radius, with 1.28 x 10 5 
particles. There was less than 3 percent difference in the orbital energy change between the lower and higher resolution 
simulations, with the lower resolution yielding a less tightly bound binary. Since the lower resolution displayed less 
dissipation, the simulations performed here can be viewed as slightly conservative in the capture rate estimations. If 
we were to allow fragmentation in the disk, a significantly higher number of particles would be required to capture the 
details. Because the energy transfer in these simulations is dominated by the large scale behavior of spiral arms and 
trailing shocks, the relatively low resolution used is sufficient. 

The impactor is initially at a distance of 1000 AU, on an orbit corresponding to a relative velocity at infinity of 2 
km s _1 . The simulations all have the same relative velocity, chosen to be representative of typical velocity dispersions 
in young clusters. The eccentricities of the encounters that result in a bound system are all less than ~ 1.1; above 
this the velocities are too high for dissipation via disk interactions to create a bound system. Tests suggest that the 
change in orbital energy is not a strong function of eccentricity for such low relative velocities. We ran prograde 
and retrograde encounters with a 1 Mq impactor using relative velocities of 0.1, 2, and 4 km s _1 , corresponding to 
eccentricities of 1.000, 1.058, and 1.235 respectively. For the prograde case, the change in energy of the different cases 
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Fig. 1. — An example of an encounter, with m = 3.0 Mq, r p = 175 AU, and i = 45 degrees. Plotted is the logarithm of the gas column 
density. Left: The impactor has not yet passed through the disk, which is in its relaxed, initial state. Right: The passage through the disk 
is completed, and spiral features are visible. 



agree to within 5 percent. The retrograde case exhibits a slightly larger variance, approximately 10 percent, with the 
highest eccentricity encounter differing from the others. By simulating a modest but non-zero velocity at infinity, we 
perhaps underestimate the capture rates for the low velocity tail of the distribution, and overestimate the high velocity 
capture rate. Since the dominant contribution to the capture rate comes from relative velocities < 2 km s _1 for all 
masses, the net contribution of this simplification is not thought to significantly alter the capture rate estimations. 

We consider four different masses m for the impactor: 0.3, 1.0, 3.0, and 9.0 M Q . We cover the range of inclination 
angles [0,180] in 45 degree increments, and we simulate encounters with unperturbed periastra r p of 0.1, 0.35, 0.6, 0.85, 
and 1.1 r<j. This provides a grid of 100 encounters. At the end of the simulation we sort the gas particles into three 
groups: those bound energetically to the impactor, those ejected from the system, and those bound to the primary. 
The stars and their associated gas particles are treated as a point mass at the center of mass of the grouping, and these 
masses are used to determine an orbital energy. A representative encounter is shown in figure ^ with m = 3.0 Mq, 
r p = 0.35rd, and inclination angle i = 45 degrees. In the right panel, the impactor has finished its passage through 
the disk. We run each encounter until the orbital energy has stopped evolving. The results of these simulations are 
summarized in figure 

The change i n energy for a give n inclination gener ally decreases with increasing periastron distance, a result in 
agreement with H ellerl l)1995|) and iBoffin et alJ (0-998). There are also differences between the low mass and high 
mass impactors, best seen by comparing the change in energy associated with prograde versus retrograde collisions for 
r p = O.lr^. Note that the disks con sidered h e re are gaseous; several authors have considered disk interactions with 
purely gravitational effects included. iPfalznerl l)2003j) showed that for low mass disks and non-penetrating en counters, 
ignori n g hydrodynamica l forces is a valid simplification, and uses this to explore a large parameter space in IPfalznerl 
(2004). Hall et al. (1996) used a reduced three-body approach to explore encounters in the absence of hydrodynamics 
or disk self gravitation, to identify the regions of the disk that dominate energy transfer. 

However, with more massive disks and penetrating encounters, orbital energy loss i n retrograde encounters is dom- 
inated by the dissipat ion of energy in the trailing shock behind the intruding star (jHellerl 119951: IBoffin et al.lll998l: 
iMoeckel fc Ballvll2006|) . As the mass of the impactor increases, this mechanism becomes more effective. For m = 0.3 
Mq, the prograde encounter is over three times more dissipative than the retrograde; for m = 9.0 M , the retrograde 
case is over twice as dissipative as a prograde passage. As the periastron radius increases, this mechanism plays less 
of a role, until at r p = Llr^ it is less effective than the spiral arms set up in a prograde encounter. The capture of 
a mass by the impactor plays a small role in the change of orbital energy. The lower mass impactors, m < 1.0 M , 
capture very little mass from the disk. The more massive impactors can increase their mass by as much as 10 percent 
during a passage, but the change in orbital energy is typically an order of magnitude larger than can be accounted for 
by simply adding mass to the impactor in a momentum conserving fashion. 
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Fig. 2. — Change in energy as a function of inclination for each periastron radius r p and impactor mass m. The diamonds are the 
simulation data points, and the curves are the interpolated values used in calculating binary formation rates. The energy is in code units, 
where mass is in M©, velocity is in km s — 1 , and distance is in AU. 



3. DETERMINING CAPTURE RATES 

In order to calculate the efficie ncy of disk-ca pture as a binary formation mechanism, we largely follow the method of 
iHellerl l|1995j) . also used bv lBoffin et alJ l|1998j) . We modify the method to account for the difference in mass between 
the impactor and the primary, and allow for a mass dependent velocity dispersion. We describe the method here. 

Determining binary ca pture rates is closely related to the calculation of stellar collision rates found in 
IBinnev fc Tremaind l|1987j) . Instead of determining the frequency of encounters with closest approach less than a 
stellar radius, we are interested in the frequency of passages that have a closest approach within some critical radius 
r c , within which the energy dissipated via disk interactions is sufficient to form a binary. This radius depends on the 
stellar masses, initial orbital energy, and inclination angle of the encounter. For a cluster of equal mass stars in a 
Maxwellian velocity distribution, the rate of encounters with a closest approach radius less than r is given by 

j BT = ^1 + , (5) 

with n the number densit y of the cluster, to* the mass of the stars and a the one dimensional velocity dispersion 
ijBinnev fc Tremainelll987ft . 

The cross section for such an encounter is easily generalized to the case of uneven mass stars M and to, with velocities 



velocities V and v, 



E(r,M,v re l) = Trr 2 + 



2-kGMt 



J rel 



(6) 



with M. = M + m the total mass of the system, and v re i = |V — v|. For calculating capture rates, the relevant 
encounter radius is not fixed, but depends on to, v re i, and inclination angle i. Via interpolation of the data in figure 
El we can find a critical radius r c , such that the encounter satisfies 



with n the reduced mass of the system, Mm/ AA. The cross section for capture is then defined as 

y(i,M,v re i) = T,(r c ,M,v re i). 

The capture rate for a star with mass M and velocity V is then obtained by multiplying the cross section by v re i 
and integrating over the velocity distribution of the encounter partners /(v); an average over the primary velocity V 
and inclination angle i is obtained by 



(7) 



(8) 



r c (m) 



Jo°° /o°° £ f(V)f(vMvrei,i,M)vreisin(i)di d 3 V d 3 V 
Jo" So /(V)«n(i)d* d?V 



(9) 



We assume a Maxwellian velocity distribution with a one dimensional, mass dependent velocity dispersion of the 
form 

where m and a are chosen to match some observations, for instance TOn = 0-3 Mm ; and an ~ 2.3 km s" 1 would be 
appropriate for the Orion Trapezium cluster ijvan Altena et alii 1 9881 Pones fc Walker! fl988|) . 

After applying this velocity distribution to equation we introduce new integration variables V = V — v, and 
v = (M a V + m a v)/(M a + to"). With this change we have 



T c (to) 



n(Mm) 



3a/2 roc pco piv 



l u Jo Jo 



exp 



(M a + m a ) 2 v 2 + (Mm) a V 2 



2to,q <7q (M a + m a ) 
y(V,i,M)Vsin(i)di d 3 vd 3 V. 



(11) 



The integral over v is straight forward, and we arrive at a formula for the capture rate of a star with mass M 
interacting with stars of mass m, 



r c (m) 



n(Mm) 3a / 2 



2nm 3 a/2 a^(M a + m a f a / 2 Jo Jo 



OO p 77 



exp 



(Mm) a V 2 



2toq (Jq (Af a + m a ) 
y(V,i,M)V 3 sin(i)di dV. 
For the case of no mass dependence in the velocity dispersion, a — 0, this reduces to 



r c (m) 



Jo 



OO /"7T 



exp 



'4a 2 



^{V,i,M)V 3 sin(i)di dV, 



(13) 
(14) 

(15) 

in equation |51 the 



which is equivalent to the formula used bv lBoffin et alJ l)1998fl . 

If instead of using the experimentally found ^(i, A4, v re i) we use the analytical S(r c , A4,v r 
inclination dependence disappears and the integral over V is easily done, resulting in a generalized collision rate, 

GM ( M 



27r(Af"+TO") 2 Q/2 

7 = 2\ : nr m n an 

' V {Mm) a ° 



1 



a^r \m 



M c 



(16) 



For the special cases a — 1 (energy equipartition) and a — (no mass dependence in the velocity dispersion), 
equation 1161 reduces to 

r2 ( 1+ GMjA 



/27rm 
7a=l = 2W — - — na r 



7a=o = 4V?i : «o'of* 2 ( 1 + 



m a 2 r J 



GM 



2a 2 r) ' 

Eauation ll8l is equivalent to equation with the mass replaced by the average of M and to. 



(17) 
(18) 



6 



TABLE 1 

The effect of the constant (3 on the 
IMF. 







P 




Mass Range 


0.0 


0.5 0.75 


1.0 


(M ) 




(% of stars) 




0.3-0.5 


42.9 


30.1 23.5 


17.2 


0.5-1.0 


35.0 


32.5 29.2 


24.6 


1.0-9.0 


21.4 


34.3 41.5 


48.1 


9.0-20.0 


0.77 


3.08 5.77 


10.1 


m/M Q 


0.32 


0.88 1.66 


3.37 



Note. — See equation l20l for the IMF 
and definition of j3. The average masses 
rh = f™ u m^(m)dm/ f™ u £,{m)dm are in- 
tegrated over the range 0.3 - 20.0 Mq. 



Equation II II is the capture rate for a single star of mass M awash in a sea of stars with uniform mass m. To arrive 
at a total capture rate for a specific star with mass M in a cluster with density n and some mass function £(m), we 
average ea uation 1111 over the mass range m; to m u , 

_ C," t{m)T c (m)dm 
J™, Z{rn)dm 

In this work we use the IMF from lKroupal pOOlj) . 

( Ci m-°- 3+l3 , 0.01 < m/M < 0.08, 
i{m) = I C 2 m- 13+/3 , 0.08 < m/M Q < 0.5, (20) 
[ C 3 m- 2 - 35 +^, 0.5 < m/M Q < M max , 

where we have added a constant (3 to all the exponents; this constant is included to approximate the effects of mass 
segregation in a cluster. The constants Ci enforce continuity at the transitions between the different power laws, and 
are functions of (3. If (3 = 0, we have the unmodified, present day IMF. As (3 increases, skewing the mass function 
toward higher values, the average stellar mass increases and encounter rates are altered. The effect of (3 on the mass 
function is summarized in table ^ Note that this single parameter alt eration of the mass function is insufficient to 
reproduce the mass function generated by dynamical mass segregation ( Gii rkan et al.ll2004|) . We use this to illustrate 
the character of mass segregation effects in a simple fa shion. The maximum mas s M max depends on the cluster mass, 
which sets the maximum stellar mass in the cluster ijWeidner fc Kroupal 12006). While we consider the primary to 
be the most massive star in the cluster, because we only simulated encounters with companions up to 9 M Q we set 
M max = 9 Mq. We now have the binary formation rate averaged over the IMF, and as a function of the level of mass 
segregation and mass dependent velocity dispersion. 

Figure [3] shows the capture rates, with no mass dependence in the velocity dispersion, for the cases where all the 
cluster members (besides the capturing star) are 0.3, 1.0, 3.0, and 9.0 M Q , as well as the rate averaged over the IMF. 
Because of the range of masses we simulated, we take mi = 0.3 M and m u = 9.0 M Q . At the high end of the mass 
spectrum this truncated mass function is fairly complete for (3 < 0.5, as stars with masses m > 9.0 M Q contribute < 3 
percent of the total population. For higher values of (3 we are missing a significant portion of the mass function, which 
leads to an underestimation of the captur e rates. The rates are normalized t o a clu s ter w ith stellar density 114 = 10 4 
pc~ 3 , and can be compared to figure 8 in lBoffin et alJ l|1998|) and figure 5 in IHellerl l|1995|) . The IMF integrated rate 
calculated here is roughly an order of magnitude higher than those previously calculated for encounters between lower 
mass stars, with mass ratios M/m ^ 0.5. 

By ignor ing any encounters which take place outside the range 0.1 r c i < r p < 1.1 we underestimate the rates 
somewhat: boffin et alJ ()1998ft find that some encounters out to r p = 2.0 result in a binary. These weak encounters 
are most important at low relative velocities. Since we are most interested in scenarios with <tq ~ 2.0 km s _1 , the 
exclusion of encounters with periastra greater than those simulated is neg l igible. The frequency of encounters with 
r p < 0.1 rd is small enough that ignoring them is a small effect. IHellerl l|1995|) finds that ext rapolating to direct 
collisions from r p = 0.2 results in changes to the capture rate on the order of 10 percent, and iBoffin et alJ l)1998fl 
reports similar results extrapolating inward from r p — 0.5 r^. 

4. DISCUSSION OF RESULTS 

4.1. Comparison to Previous Work 

In this section we consider the case of no mass dependence in the velocity dispersion, and no mass segregation 
in the IMF, i.e. a = f3 = 0. Table [5] shows the capture rates from this work, for a variety of star forming en- 
vironments, compared to rates calculated in previous work. For the Trapezium center we take n = 2 x 10 4 pc -3 
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Fig. 3. — Binary formation rates normalized to a cluster with number density m = 10 4 pc~ 3 . Shown are the rates for a 2 M q star 
with a 2 Mq, 500 AU disk in a cluster of stars with uniform mass 0.3, 1.0, 3.0, and 9.0 Mq, and in a cluster the IMF of eguation l20l for the 
range 0.3 Sm/M® ^ 9.0. In this plot there is no mass dependence in the velocity dispersion, nor mass segregation in the mass function. 



TABLE 2 

Capture rates in different environments. 



Region 


n 

(PC" 3 ) 


a 

(km s" 1 ) 


This Work 


References 
1 

r c 

(Myr) 


2 


Trapezium, center 


2 x 10 4 


2.3 


0.53 


0.06 


0.018 


Trapezium 


2 x 10 3 


1.5 


0.17 


0.02 


5 x 10~ 3 


Open cluster 


10 2 


1 


0.02 


3 x 10" 3 


8 x 10" 4 



References. — (1) Boffin et al. 1998; (2) Heller 1995. 



l|Hillenbrand fc Hartmannlll 998ft. and g n = 2.3 km s" 1 i|van Altena et al]ll988H.Tones fc Walked! 1988ft . The other re- 
gions considered are those in lBoffin et all (fT998). For a massive star in a cluster of lower mass companions, the capture 
rates are roughly an order of magnitude higher than encounters involving impactors with masses M/2 < m < 2M. 

Figure 0] shows the probability of an average encounter resulting in a captured binary for each of the impactor 
masses. The probability is given by T c (m)/T enc , with T enc the total encounter rate for r < 1.1 found by averaging 




equat ionll8|o ver th e IMF. The probabilities are one to two orders of magnitude larger than those of lHellerl l)1995j) or 
Boffi n et al.l (^998). However, despite the lowest mass impactors, m, = 0.3 M , encountering a disk nearly an order 
of magnitude more massive than themselves, they have the lowest probability of capture. The increased gravitational 
radius of the more massive impactors is enough to overwhelm any effect such a relatively massive disk might have on 
a low mass impactor. 

The increased effectiveness of disk interactions in capturing companions seen here is due to two different effects. 
Due to their greater gravitational focusing, massive stars in a cluster of lower mass companions will have a larger 
encounter rate, which though not the main reason for the larger rates here, contributes somewhat. The main cause 
of these higher rates lies in the ratio of the disk mass to its radius. The primary-disk systems in this work have the 
same ratio md/m p = 0.1 as Heller (1995). However, while these disks are 20 times as massive as Heller's, the radii are 
only 4 times as large. The disks in Boffin ct al. are both less massive and have larger radii. It is this higher relative 
concentration of disk material that is the dominant cause of the high rates calculated here. 

In lower density regions such as an open cluster, the increased efficiency of high mass capture is insufficient to make it 
a significant contributor to the total binary fraction. However, in dense regions such as the Trapezium cluster in Orion, 
with stellar density ~ 2 x 10 4 , one could expect more than 50% of the massive stars to have captured a companion 
after 1 Myr. This is a high enough rate to be considered a real contribution to the observed multiplicity of high mass 
stars. 

In the Arches cluster near the Galactic center, the cen tral mass density is ~ 3 x lOfMg pc" 3 llFis-er et al ] ll9991. 
compared to the Trapezium core at ~ 3 x 10 4 M© pc -3 (Hil lenbrand fc HartmannH T998). At such high densities the 
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capt ure rate T r could be well above unity after 1 Myr, though the only limit on the velocity dispersion is a < 22 km 
s" 1 l)Figer et all 120 02). which would drive the rate down to negligible values. It is worth noting that the majority of 
encounters that result in capture occur at low relative velocities v re i ^ 4 km s . As the velocity dispersion increases 
and such encounters are less frequent, the capture rate plummets, clearly seen in figure 

4.2. Effect of a Mass Dependent Velocity Dispersion 

Shown in figure[5]are the IMF integrated capture rates, for different values of a in eauation llOl as well as the velocity 
dispersions as a function of mass. The results below depend not only on the value of a but on mo, the anchor point 
of the velocity dispersion. We choose m = 1.0 M Q so that the dispersion on the low end, rn* ~ 0.3 M , is not 
unrealistically high, yet a good differentiation between the most massive and least massive stars' dispersions can be 
seen. 

As the cluster tends toward equipartition, a = 1, the binary formation rate roughly triples. The cause of this is 
that the massive stars are moving more slowly, and the lower relative velocity increases the encounter rate and the 
probability of capture. For the case of a = —0.25, with the more massive stars moving faster, the rate drops by roughly 
a factor of two. 

A cluster with stars in energy equipartition will tend to have more captured binaries involving its most massive stars. 
However, in most clusters equipartition is by no means a good assumption. Observationally, ther e is evidence in Orion's 
Trapezium cluster for a weak tr end toward increasing velocity dispersion with decreasing mass Ipones &: Walkerlll988t 
Hillenbrand &: Hartmaimlll998|) . Hillenbrand & Hartmann find a = 2.81 km s" 1 for 0.1 < M/M Q < 0.3, and a = 2.24 
km s" 1 for 1.0 < M/Mq < 3.0. In our formula for the mass dependence of the velocity dispersion, a value of a ~ 0.2 
seems appropriate for the Trapezium. At such low values, the increase in the capture rate is only ~ 25 percent, much 
less dramatic than a cluster in equipartition. 

4.3. Effect of Mass Segregation 

An effect which would raise the rates c alculated here is m ass segregation; there is evidence for this in Arches and 
Trapezium clusters. In the Arches cluster. TStolte et"al~l (|2002D find that the overall slope of the IMF flattens inside 5", 
from a value of — 2.7 ± 0.7 to —1. 1 ± 0.3. In Trapez i um, the mass function of the clu ster taken as a whole (r < 2.5 
pc) is similar to a standard IMF f Hille nbrand! 1 1 997}) . Hillenbrand fc Hartmanr] l)1998|) find strong evidence for mass 
segregation toward the cluster core for masses down to ~ 1 M©, and find that within 0.05 pc of the cluster core the 
mean mass per star fh ~ 2 — 4 M©, compared tom~ .75 M© for r > 0.1 pc. We approximate the effects of segregation 
by increasing the value of f3 in equation[5D| In order to raise the average stellar mass above 1 M©, we need (3 ~ 0.75, 
though for the Arches cluster /3 = 1 would be more appropriate. 

As (3 increases from 0, the more massive stars become over-represented in the mass function compared to the present- 
day IMF. In the absence of a mass dependent velocity dispersion, the total binary formation rate essentially follows 
the rate that would be found if all the stars had the average stellar mass. What is more interesting is the capture rate 
broken down by mass. Shown in figure H3 are the mass functions for different values of f3, and the capture rates as a 
function of mass weighted by the mass function. The velocity dispersion ctq = 2 km s -1 . 

As (3 approaches increases from to 1, the mass function flattens and 9 M© stars increase their representation 
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in the stellar population by over an order of magnitude. For (3 — the ratio £(0.5)/£(9.0) ~ 850. For 0=1, 
£(0.5)/£(9.0) ~ 50. This change in the proportion of massive stars is reflected directly in the capture rates. Since 
massive stars have higher capture probabilities (see figure the ratios of capture rates are smaller. For (3 = the 
rate for 0.5 M impactors is ~ 50 times as large as the rate for 9.0 M Q stars; at j3 = 1 the rate for 0.5 M stars is 
only ~ 7 times as high as that of 9.0 Mq stars. If mass segregation occurs in a young cluster on a timescale short 
enough that disk capture has a chance to contribute to the multiplicity of massive stars, presumably found near the 
center of the cluster, than capturing relatively massive companions would be likely. In an environment like the Arches 
cluster, where the mass function is in fact flat, a massive companion is the most likely outcome of a capture event. 
Mass segregation thus has the double effect of increasing the net capture rate, and biasing the mass of captured stars 
toward the high end of the mass spectrum. 

4.4. The Fate of the Binaries 

With semi-major axes of several x 10 3 AU, the resultant binaries from these captures present extremely large cross 
sections for further interactions with other cluster members, especially in the dense (n ~ 10 3 — 10 4 pc -3 ) environment 
of a young, rich star cluster. If binaries formed via disk encounters are to contribute to the observed high multiplicity 
of massive stars, they need to be robust enough to survive. 

A full treatment of this question would include a campaign of n-body simulations, beyond the scope of this paper. In 
order to get some sense of the fate of captured binaries, we instead compare the energy of an encounter with another 
cluster member to the orbital energy of the binaries. In considering an encounter between a binary and a single star 
there is a critical rela tive velocity, above which the system's total energy is positive, and total ionization is possible 
(see e.g. lHeggieiri975l) . For a single star of mass mi encountering a binary with total mass Ai, reduced mass \x and 
semi-major axis a, this relevant comparison of energies is 

1 Mm 2 2 _ GMn 

2 M + m 2 Vc ~ 2a ' [ ' 

The term on the left is the energy associated with the encounter between the field star and the binary, moving with 
relative velocity v c . The term on the right is the negative of the orbital energy of the binary. For relative velocities 
greater than v c , ionization is energetically permitted. 

For velocities below v c , these encounters can be resonant, with complex orbit s eventually leading the ejection of 
one of the stars, or a rela tively simple exchange wherein a new binary is formed ijHut fc Bahcal]|IT983t iHeggie fc Hutl 
2003: iFregeau et al"ll2004|) . The critical velocity v c is not a hard boundary, but for velocities above ~ 2v c , ionization 
is the dominant result. In the case of an encounter which forms a temporary, unstable three body s ystem, the least 
massive member is preferentially ejected l|Anosoval 198a iMikkola fc V altoncn 19 86llBate et al.ll2002|) . Note, though, 
that this "dynamical biasing" is not a rule, especially with soft binaries l)Heggidll975t) . Additionally, most numerical 
studies of binary-single star interactions focus on nearly equal mass components, and it is not immediately clear that 
the results apply to a situation with one dominant mass. However, as a rough guide we continue the analysis, with 
the assumption that some (if not most) of the time, an encounter with another cluster member will leave the most 
massive star in a binary. 

Shown in figure \7\ are these energies for various combinations of companion mass m and the mass of the second 
collision partner m%. The lines are the left hand term in equation 1211 for different values of m.2, with the relative 
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Fig. 7. — Energies associated with the final binary orbits (points), and encounters with other cluster members of various masses 7712 
(lines). Plot titles indicate the mass of the binary companion to the 22 Mq primary. The absolute value of the negative binary e nergy is 
plotted for comparison with the encounter energies. Points above the solid lines are energetically unable to be ionized. See section [4.41 for 
further details. 



velocity equal to the most likely relative velocity, twice the velocity dispersion. The points are the results from a 
Monte Carlo sampling of the simulation results. For a random velocity dispersion, the relative velocity is drawn 
from a Maxwellian distribution with a = \flo§. The inclination angle i and periastron r p are sampled subject to 
their probabilities, proportional to sin(i) and r p respectively. The sampling of relative velocity, inclination angle, and 
periastron is repeated until an encounter results in a binary. In figured the condition for a binary to be more likely 
to survive is that its position is above the curves, i.e. the absolute value of the binary's orbital energy is greater than 
the energy of the encounter with a third body. 

Capture formed binaries with low mass companions, m ~ 0.3 Mq, appear unlikely to survive for long in a dense 
environment. Even encounters with slow moving, low mass cluster members are likely to ionize a typical binary. As the 
mass of the captured companion increases, more of the points lie above the energy curves. For m > 3.0 Mq, encounters 
with low mass cluster members at moderate velocities are largely unable to ionize the system, but the majority of 
binaries seem soft. An effect not taken into account here is further passages through the disk. If the orbital period of 
the binary is less than the encounter timescale in the cluster, the binary will harden before encounters with further 
cluster members. In practice this effects only the more massive binaries; plotted in figure [3] are the orbital periods 
versus encounter times of the points in figure for m — 9.0 Mq. The encounter times were calculated from equation 
1181 using the average mass of a cluster member ~ 0.5M Q , and the semi-major axis of the binary as the encounter 
radius. If the orbital period is short compared to the encounter time, the point will lie above the solid line. The 
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Fig. 8. — Encounter times versus orbital times for the same points as figure [7| for m = 9.0 Mq. Points above the lines will encounter 
the disk a second time before the binary encounters another star in a cluster with the labeled stellar densities. 



encounter time was calculated for n — 10 4 pc~ 3 , and adjustments for different stellar densities are plotted. 

To estimate the effect of repeated disk passages, if a point lies above the line we assume that it passes through 
the disk a secon d time. The change in energy associated with further passages is dependent on inclination angle and 
impactor mass l|Moeckel fc B allv 200(|); as an approximation we assume that the second passage results in one half 
the energy change as the first encounter. The survival percentages of the binaries are shown in figure The figure is 
the result of Monte Carlo sampling of the binaries that form, and their encounters with further cluster members. At 
each velocity dispersion, 4000 binaries are created, with initial conditions chosen as detailed above. After each binary 
is formed, an encounter partner is drawn from the same velocity distribution and from the IMF of equation 1201 The 
energy of the binary is compared to the energy associated with the encounter. If the encounter energy is insufficient to 
ionize the binary, we count the binary as a survivor. Plotted is the number of survivors, Ns, divided by the number of 
binaries tested, Nb- Poisson counting statistics are assumed for the error bars. Plotted for each mass are the results 
with and without accounting for second passages, in a stellar environment similar to the Trapezium center, n = 2 x 10 4 
pc~ 3 . 

When further passages through the disk are included in the calculation, there is a modest increase in the survival of 
high mass binaries, and essentially no change to the lower mass companions' survival rate. The chances of a low mass 
(m < 1.0 Mq) companion surviving in a binary after further encounters with cluster siblings are low for all but the 
lowest velocity dispersion environments. In a region like the Trapezium, with oq 2.3 km s _1 , the survival percentage 
for 1 Mq companions is ~ 0.05. However more massive companions, with m > 3.0 M Q , have a much higher survival 
percentage. A 9 Mq companion cannot be ionized by ~ 60 percent of further encounters in the same environment. 
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Fig. 9. — The number of surviving binaries Ng after an encounter with another cluster star, divided by the number of binaries tested, 
Nb, based on Monte Carlo sampling of the cluster. Solid lines show the survival percentage without accounting for second encounters, 
while the dotted lines show the percentage with second encounters included for n = 2 X 10 pc -3 . The curves are shown for the labeled 
companion masses m. 



It seems that if stellar encounters are frequent enough for disk assisted capture to be a significant binary formation 
mechanism, further encounters will tend to disrupt the lower mass binaries, with higher mass binaries preferentially 
surviving. An effect that this analysis does not explore is the progressive softening of a binary by repeated, non-ionizing 
encounters. If an interaction results in a less tightly bound binary, the next encounter will be more likely to disrupt 
the system. This is beyond the simple energetic treatment given here, but has the potential to further decrease the 
survival rates of binaries near the critical binding energy. 

4.5. Observable Indications 

Binaries formed by capture would exhibit several observable characteristics. The semi-major axis a of the orbit 
is given by a = GAifi/2E, with E the orbital energy. Typical values for the encounters simulated here range from 
a ~ 10 4 AU for m — 0.3Mq to a ~ 10 3 AU for m — 9.0Mq. Since the periastra of the captured binaries are less than 
the disk radius, the eccentricities are large. Thus the capture of an initially unbound impactor results in a long period, 
highly eccentric (e ~ 0.9) binary. Previous w ork on repeated encounters between a massive star-disk system and a 
less massive impactor l)Moeckel fc Ballvll2006j) suggests that the disk will be truncated within a few orbits to roughly 
the periastron distance, as the semi-major axis of the binary shrinks. During these repeated encounters the orbit is 
circularized to some extent, though it remains eccentric with 0.6 < e < 0.8 depending on the inclination angle. The 
evolution of the orbit's semi-major axis is likewise inclination dependent, shrinking by anywhere from a factor of ~ 3 
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for a prograde encounter to an order of magnitude for retrograde. In the most dissipative cases, in-plane collisions, 
the disk is almost totally disrupted. In out-of-plane encounters, the disruption of the disk is more gradual, with the 
disk retaining material through ~ 10 orbits or more. 

The disruption of the circumstellar environment during capture could help explain the observed lack of disks around 
massive stars in the Trapezium core and other young clusters. As opposed to young, embedded protostars, where there 
are examples of disks, there is no evidence of circumstellar disks around mature massive stars. In the Trapezium, the 
most massive stars are multiple, and diskless. A higher mass binary companion is more spatially disruptive to the disk 
during repeated passages, and is also more likely to have time to pass through several times before an encounter with 
another cluster star. Regardless if the ultimate fate of the binary is survival or disruption, the initial passages through 
the disk provide an efficient disk disruption mechanism. Since binary capture is destructive to the disk, the disk 
dissipation timescale is tied to the binary formation timescale. In dense clusters, n > 4 x 10 4 pc~ 3 , the timescale for 
all massive stars to have undergone a binary capture is under 10 6 yr, comparable to the length of the embedded phase, 
suggesting that destr uction via repeated en counters could contribute to the dissipation of the disks. This is similar 
to an effect noted bv iPfalzner et all l)2006j) . though in their scenario disk destruction around massive stars occurs 
via many encounters with different low mass stars, whereas here it is due to repeated encounters with a captured, 
moderate mass companion. 

The tilting and warping of a disk in an out-of-plane encount er has been seen in many previous simulations l)Helleil 
119951 lLarwo^ll997t iBoffin et alJll99St iMoeckel k Ballvl l2"006) . Shown in figure ITU1 is the change in orientation, with 
respect to a fixed coordinate system, of the disk due to the passage of the impactors. The same mass, inclination, 
and periastra as in Figure [3 are plotted. The disk orientation is calculated by averaging the angular momenta vectors 
of the particle orbits about the primary. The change is given in degrees, and is simply the angle between the initial 
and post-passage orientation vectors, not precession about the axis of the binary orbit. The disk is warped by these 
passages, so the tilt of the disk varies with radius; plotted here are the orientations measured by considering disk 
material within 50 AU of the primary. 

The change in disk orientation from an encounter with an impactor of mass m < 1.0 M is not likely to be more 
than a degree. However, for the larger impactor masses simulated, changes in orientation of ~ 5 degrees or more are 
possible, with the most extreme encounters twisting the disk through 10s of degrees. During out-of-plane encounters, 
the remnant disk experiences an almost impulsive torque during periastron passages. Assuming that the orientation 
of an outflow is somewh at tied to the angular momentum of th e disk that feeds it, such periodic changes could be 
indications of a capture l)Bate et al.ll200ft |Moeckel "fc"B allv 20Q^). 

The statistics of mass ratios q = m/M in high mass multiple systems are not constrained to the point where 
a comparison to this formation mechanism can be made. Because higher mass companions preferentially survive, 
one would expect the distribution of mass ratios to be biased toward the higher end. For visual binaries in Orion, 
iPreibisch et al.l l)1999[) concluded that the di stribution of q was consistent with the companions being drawn from 
the IM F. In NGC6611. iDuchene et all ((2001) could neither confirm nor rule out an IMF distribution. iMason et alJ 
(1998) surveyed spectroscopic binaries among O stars, finding a flatter q distribution and a lack of low q systems, 
however these results may be partially due to bias against detecting low mass companions in spectroscopic searches. 
The tendency of the binaries studied here to be relatively high q seems to be currently unconstrained by observations. 

5. SUMMARY 

The capture of companions through disk interactions has been previously studied in the context of low mass stars. 
In this paper we have extended the examined parameter space toward values more representative of high mass stars. 
The capture rates we calculate are an order of magnitude or more higher than those found in the lower mass scenario. 
Given the current stellar density and velocity parameters of a region like the Trapezium core, the rates calculated here 
could account for ~ 50% of the massive stars having a companion after 1 Myr. Though not high enough to account 
for the extremely high companion fraction of the Trapezium stars, this rate is not negligible, and the disk-capture 
mechanism should not be ignored in dense regions of massive star formation. However, the binaries that form are soft 
and perhaps easily disrupted. 

We attempted to estimate the robustness of the binaries, concluding that massive companions are more likely to 
survive in a dense cluster. Lower mass companions tend to form very weakly bound binaries, which are unlikely to 
remain bound after encounters with other cluster members. The same conditions that lead to a high capture rate, 
low velocity dispersion and high stellar density, result in a high rate of disrupting encounters. If this mechanism is 
to contribute to the observed multiplicity, the most likely surviving binaries will have massive companions, and will 
have passed through the disk several times before another encounter occurs. However, the details of binary-single star 
encounters, with one dominant mass in an eccentric binary and the remnants of a disk still around, warrant further 
investigation beyond the simple treatment given here. 

The effects of a mass dependent velocity dispersion are not likely to affect the capture rates, except for extreme cases, 
close to true energy equipartition. Mass segregation, however, enhances the capture rates significantly for more massive 
impactors. In a heavily mass segregated environment, such as the Arches cluster, capture of a massive companion 
(m > 5 Mq) would be a frequent occurrence, though in the Arches cluster itself the high velocity dispersion would 
lower the capture rate significantly. Capturing a companion is destructive to the disk, and the timescale for binary 
capture is comparable to the expected disk lifetime in dense regions. The likelihood of a companion being massive, 
and the disk destruction associated with repeated encounters, could help explain the observed lack of disks in regions 
like the Trapezium. 
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This work has ignored the effects of disk replenishment from the reservoir of gas at the center of a young cluster, 
and the dissipation of the disks over time. Dissipation can occur viscously or from radiative mechanisms, and the 
replenishment of the disk depends sensitively on the details of the radiation field of the massive star. Whether these 
combined effects enhance or decrease the rates calculated here will depend on the age of the protostars, the local gas 
and dust structure, and the amount of gas left in the cluster. 
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